Magnetic and thermodynamic properties of Li2VOSi04 : a 
two-dimensional S =1/2 frustrated antiferromagnet on a square 

lattice 



R. Melzi, S. Aldrovandi, F. Tedoldi and P. Carretta* 
Dipartimento di Fisica "A. Volta" e Unitd INFM di Pavia, Via Bassi 6, 27100 Pavia, Italy 

P. Millet 

Centre d'Elaboration des Materiaux et d'Etudes Structurales, CNRS, 31055 Toulouse Cedex, 

France 

F. Mila 

Institut de Physique Theorique, Universite de Lausanne, CH-1015 Lausanne, Switzerland 

(February 1, 2008) 



Abstract 



NMR, yuSR, magnetization and specific heat measurements in Li2VOSi04 
powders and single crystals are reported. Specific heat and magnetization 
measurements evidence that Li2VOSi04 is a frustrated two-dimensional S = 
1/2 Heisenberg antiferromagnet on a square lattice with a superexchange 
coupling Ji, along the sides of the square, almost equal to J2, the one along 
the diagonal (J2/J1 = 1.1 ± 0.1 with J2 + Ji = 8.2 ± 1 K). At ~ 2.8 
K a phase transition to a low temperature collinear order is observed. Tc 
and the sublattice magnetization, derived from NMR and /iSR, were found 
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practically independent on the magnetic field intensity up to 9 Tesla. The 

critical exponent of the sublattice magnetization was estimated /3 2± 0.235, 

nearly coincident with the one predicted for a two-dimensional XY system 

on a finite size. The different magnetic properties found above and below Tc 

are associated with the modifications in the spin hamiltonian arising from a 

structural distortion occurring just above T^. 
PACS numbers: 76.60.Es, 75.40.Gb, 75.10.Jm 
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I. INTRODUCTION 



In recent years one has witnessed an extensive investigation of quantum phase transition 
in low- dimensional S = 1/2 Heisenberg antiferromagnets (QHAF) as a function of doping, 
magnetic field and disorderB. For example, two-dimensional QHAF (2DQHAF) have been 
widely studied in order to evidence a phase transition from the renormalized classical to 
the quantum disordered regime upon charge dopingi. Another possibility to drive quantum 
phase transitions in a 2DQHAF is to induce a sizeable frustration. In particular, for a square 
lattice with an exchange coupling along the diagonal J2 about half of the one along the sides 
of the square Ji (see Fig. la), a crossover to a spin-liquid ground state is expectediH^. 
For J2I J\ ~ 0.35 Neel order is realized, while for J2I J\ ~ 0.65 a collinear order should 
develop. The collinear order (see Fig. la), which can be considered as formed by two 
interpenetrating Neel sublattices with staggered magnetization ni and n2, is characterized 
by an Ising order parameter a = ni.n2 = ±ll. The two values of a correspond to the 
two collinear configurations which can develop at low temperature (hereafter T), one with 
spins ferromagnetically aligned along the x axis, corresponding to a magnetic wave-vector 
Q = (0, vr/a), the other with spins ferromagnetically aligned along the y axis (Q = (vr/a, 0)). 
At a certain temperature a phase transition occurs and the system choses among the x or y 
collinear configurations. The precise boundaries of the J2/J1 phase diagram for a frustrated 
2DQHAF are unknown and could be modified by the presence of a finite third neighbour 
coupling!. Most of these theoretical predictions have not found an experimental support so 
far, mainly due to the absence of systems which can be regarded as prototypes of frustrated 
2DQHAF on a square lattice with J2 close to Ji, even if some frustrated 2DQHAF with 
different topologies have been recently studied!. Moreover, a theoretical description of the 
spin dynamics of these frustrated 2DQHAF as a function of T and magnetic field intensity 
is still missing. 

Recently, two vanadates which can be considered as prototypes of frustrated 2DQHAF on 
a square lattice have been discoveredi: Li2VOSi04 and Li2VOGe04. These two isostructural 
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compounds are characterized by a layered structure containing V"^^ (S* = 1/2) ionsi (see Fig. 
lb). The structure of V^^"*" layers suggests that the superexchange couplings between first 
and second neighbours are similar. It is however difficult a priori to decide which one 
should dominate: first neighbours are connected by two superexchange channels, but they 
are located in pyramids looking in opposite directions and are not exactly in the same plane, 
whereas second neighbours are connected by one channel, but are located in pyramids looking 
in the same directions and are in the same plane. On the basis of NMR and susceptibility it 
has been possible to demonstrate that in Li2VOSi04 J2I J\ is of the order of unity and that 
the ground state is collineari, as expected for J2I J\ ~ 0. esi. 

In this paper we present a detailed study of the magnetic and thermodynamic properties 
of Li2VOSi04 by means of NMR, /iSR, magnetization and specific heat measurements. In 
particular, we show that the spin dynamic and static properties above the collinear ordering 
temperature Tc are consistent with the ones theoretically predicted for a frustrated 2DQHAF 
with J2I J\ — 1- The phase transition to the collinear phase seems to be triggered by 
a structural distortion occuring just above T^, which possibly modifies the superexchange 
couplings and lifts the degeneracy among the two ground state configurations. The critical 
exponent of the sublattice magnetization and the independence of Tf. on the magnetic field 
intensity up to 9 Tesla, suggest that the transition is driven by the XY anisotropy. 

The paper is organized as folows: in Section II we present the experimental results 
obtained by each technique; in Section III we discuss the experimental results in the light of 
numerical and analytical results for frustrated 2DQIIAF on a square lattice, first for T > 
and then for T < T^, the main conclusions are summarized in Section IV. 

II. EXPERIMENTAL ASPECTS AND EXPERIMENTAL RESULTS 
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A. Sample Preparation, Specific Heat and Magnetization 



Li2VOSi04 was prepared by solid state reaction starting from a stoechiometric mixture 
of Li2Si03, V2O3 and V2O5 according to the procedure described in Ref. 9. The sample 
was analyzed by X-ray powder diffraction (XRD) using a Seifert C3000 diffractometer with 
CuKa radiation and then pressed into a 1 g pellet followed by a short sintering in vacuum 
at 800 C for 6 hours. Single crystals, of average size 1 x 1 x 0.2 mm^, were obtained from 
Li2VOSi04 powder heated at 1150 C for 2 hours, slowly cooled at a rate of 5 C/hour down 
to 1000 C and then furnace cooled down to room temperature. 

Specific heat (C(T)) measurements have been performed on a sintered pellet of 
Li2VOSi04 by using a standard homemade adiabatic calorimeter. The contribution of the 
addenda decreased from about 5% to below 1% of the total heat capacity on decreasing T 
from 25 K to 2.5 K. At low T the specific heat shows a broad maximum due to the correlated 
spin excitations and a sharp peak around 2.8 K (see Fig. 2a) associated with a second order 
phase transition, as can be inferred from the non-singular behaviour of the entropy around 
2.8 K. Above 20 K a rapid increase, originating from phonon excitations is observed (see 
Fig. 2a). In order to accurately estimate the magnetic contribution C""(T) (Fig. 2b) to the 
specific heat one has to subtract the phonon term C^{T), extrapolated to low T. Cp{T) was 
observed to follow a Debye law from 20 to 70 K, namely 



with Qu ~ 280 K the Debye temperature. It must be stressed that below 15 K C^iT) is 
negligible with respect to C"^{T), therefore, any incorrect extrapolation of C^{T) to low T 
will not affect the estimate of C"*(r) below 15 K. 

Magnetization measurements were carried out with a Quantum Design MPMS-XL7 
SQUID magnetometer, both on powders and on single crystals. The T dependence of the 
spin susceptibihty x — ^1 H is shown in Fig. 3. One observes a high T Curie-hke behaviour, 
a low T maximum around 5 K and a kink at Tc ~ 2.8 K, the same temperature where a peak 




(1) 
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in the specific heat is detected. The kink is better evidenced if one reports the derivative of 
the susceptibihty dx/dT (see the upper inset to Fig. 3). The susceptibihty above 15 K can 
be appropriately fitted by 

xiT) = xvv + cJiT + Q), (2) 

where 9 is the Curie-Weiss temperature, = N^g'^S^S + l)fi%/3kB {g the Lande factor 
and fiB the Bohr magneton) is the Curie constant and xvv the Van-Vleck term. The best 
fit of the data in the T range 10-300 K yields = 8.2 ± 1 K0, = 0.34 emu K/mole 
and Xvv = 4 x 10~^ emu/mole. We point out that the value of is in good quantitative 
agreement with what one would expect for an S = 1/2 paramagnet, while the absolute value 
of Xvv is consistent with a separation between the d^y ground state and the first excited t2g 
levels of the order of 0.15 eV, which is typical for V^^ in a pyramidal environmentS. Below 
Tc the T dependence of the susceptibility for magnetic fields H \\ c and if ± c is different, 
as expected in the presence of long range order. In particular, one observes that while for 
_L c the susceptibility progressively diminishes on decreasing temperature, for H \\ c it 
flattens (see the lower inset to Fig. 3), suggesting that V''"'" magnetic moments lie in the ab 
plane. 

The magnetic field dependence of Tc, derived either from the kink in the susceptibility 
or from the maximum in dx/dT, was measured from 0.1 up to 7 Tesla, where the Zeeman 
energy gfisH is greater than ksQ, and, remarkably, T^ was not observed to vary by more 
than 0.07 K, i.e. less than OMT^ (see Fig. 4). 

B. /iSR 

Zero field (ZF) /iSR measurements have been carried out on Li2VOSi04 powders at ISIS 
pulsed source, both on EMU and MUSR beamlines, using spin-polarized 29 MeV/c muons. 
The time evolution of the muon polarization is characterized by a constant background, due 
to the sample holder and cryostat walls, and by a fast decay which progressively changes from 
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exponential (see Fig. 5) to gaussian on increasing temperature, for T > Tc. Below Tc both 
oscillating and non-oscillating components are evident, the second one with an amplitude 
about half of the former one, as usually expected in magnetic powders with equivalent muon 
sites0. It must be mentioned that below Tc around 10% of the total asymmetry is missing, 
possibly due to fast precessing muons which cannot be detected at a pulsed muon source. 
Summarizing, below Tc the time evolution of the muon polarization was fitted according to 

P^it) = Aback + Aie-"* cos(75^t + 0) + A2e-^\ (3) 

where Aback is the sample holder background, Ai is the amplitude of the oscillating compo- 
nent, with 7 = 27r X 135.5 MHz/Tesla the fi~^ gyromagnetic ratio and the local field at 
the /i"*", while A2 is the ampltude of the non-oscillating component with A the longitudinal 
decay rate. Above Tc the polarization was fitted by 

P,{t) = Aback + Ae-''e-''-''/' (4) 

where the exponential term is the relaxation induced by the progressive slowing down of the 
V^^ spin fluctuations on decreasing temperature, while the gaussian term should originate 
from nuclear dipolar interaction. In particular, it is likely that localizes close to the apical 
oxygens, where it is coupled to ^Li nuclear magnetic moments. The gaussian relaxation rate 
was estimated = 0.34 ± 0.01 fis'^, for Tc < T < 4.2 K, a value typical for relaxation 
driven by nuclear dipole interactionlll. 

The T dependence of the local field at the muon and of the longitudinal relaxation 
rate A, derived after the fit of the data with Eqs. 3 and 4, are reported in Figs. 6 and 7, 
respectively. B^{T), which yields the T dependence of V^^ average magnetic moment, is 
characterized by a sharp but continuous decrease on approaching Tc, while A(T) is charac- 
terized by a divergence at Tc, as expected for a second order phase transition. 
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C. ^Li and '^^Si NMR 

^Li (7 = 3/2) NMR measurements have been carried out both on single crystals and 
powders, while ^^Si (/ = 1/2) NMR, due to the reduced sensitivity could be performed only 
in powder samples. The measurements have been performed using standard NMR pulse 
sequences. In particular, the spectra have been recorded by Fourier transform of half of 
the echo signal when the line was completely irradiated or by summing spectra recorded at 
different frequencies when it was only partially irradiated. The NMR resonance frequency 
of ^Li was observed to shift to high frequencies on decreasing temperature, with a trend 
identical to the one of the macroscopic susceptibility (Fig. 3). In fact, for ^Li NMR shift one 
can write 

7a/C(T) = ^^•^^^y^ + S (5) 

where Aj is the hyperfine coupling tensor with the j'-th V^"*" and b the chemical shift. A T 
dependent shift was observed both in the single crystals and in the powders, evidencing a 
sizeable transferred hyperfine interaction of V^"*" spins with ^Li nuclei. From the plot of the 
shift versus the susceptibility (Fig. 8) the hyperfine coupling constants and the chemical 
shift were determined. It turns out that 5 = —70 ± 30 ppm and that the hyperfine field at 
^Li is given by 

^ = E(^'^^p).4 + E ^ts.-. (6) 

where Adip is the dipolar coupling with V^"*" ions, while At = 850 Gauss is the transferred 
coupling, which is supposed to arise from the two V^"*" nearest neighbours only. On the other 
hand, ^^Si NMR resonance frequency in the powders is constant from room temperature 
down to 4 K, pointing out that the hyperfine coupling is of purely dipolar origin in this case. 

Below Tc, in the single crystals, for H \\ c, ^Li NMR spectrum splits in three lines (see 
Fig. 9): a central one with an intensity about twice of that of two equally spaced sateUites. 
The two sateUites correspond to ^Li sites with hyperfine fields of equal intensity but opposite 
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orientations, while the central line corresponds to ^Li sites where the hyperfine field cancels 
outi. The T dependence of the satellites shift is proportional to the amplitude of V^^ 
magnetic moment and, therefore, it is another method, besides ZF /iSR, to determine the 
temperature dependence of the sublattice magnetization (see Fig. 10). 

^^Si NMR powder spectrum shows a quite different behaviour at low T. Around 3 K, still 
above Tc, one observes the appearence of a shifted narrow peak (see Fig. 11). On decreasing 
T the low-frequency peak progressively disappears, while the intensity of the high frequency 
one increases. This jump in ^^Si NMR shift has to be associated either with a modification 
of the chemical shift or of the hyperfine coupling, suggesting the occurrence of a structural 
distortion just above T^. It has to be noticed that, on the contrary, no anomaly was detected 
in ^Li spectra around 3 K. Below ^^Si NMR linewidth is very close to the one above T^, 
indicating that the local field at ^^Si site is zero. 

Nuclear spin-lattice relaxation rate 1/Ti was measured by exciting the nuclear magne- 
tization either with a comb of saturating pulses or with a 180° pulse (inversion recovery 
sequence). Both for ^Li and ^^Si the recovery of nuclear magnetization towards equilibrium 
was a single exponential, indicating that for ^Li also the ±3/2 ±1/2 lines were size- 
ably irradiated during the measurements. In fact, at room temperature one can discern the 
±3/2 ^ ±1/2 lines shifted by ~ 40 kHz from the ±1/2 -1/20. The T dependence of 
^Li 1/Ti is shown in Fig. 12. One observes that 1/Ti is constant from room temperature 
down to ~ 3.2 K, then shows a peak at and rapidly decreases in the ordered phase. ^^Si 
1/Ti shows a similar T dependence below 4.2 K (Fig. 13), its absolute value, however, is 
about two orders of magnitude smaller, supporting the conclusion in favour of a hyperfine 
coupling of purely dipolar origin. 

III. DISCUSSION 
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A. Above Tc 



The T dependence of the susceptibihty and of the specific heat allows to derive informa- 
tion on the basic parameters of the electron spin hamiltonian, namely the coupling constants 
and their ratio J2I J\- For a non-frustrated S* = 1/2 2D Heisenberg AF on a square lat- 
tice the Curie- Weiss temperature O = Ji nearly coincides with the temperature where the 
susceptibility displays a maximum and one has = 0.93560. On the other hand, in 
Li2VOSi04 6 = J2 + Ji = 8.2 ± 1 K is significantly larger than = 5.35 K (see Fig. 3), as 
expected for a frustrated system. By comparing the measured ratio = 0.65 ±0.07 with 

exact diagonalization and quantum Monte Carlo (QMC) results it is possible to estimate 
J2/Jii. It turns out that J2/J1 is close either to 0.25 or 2. 5!, however, it is not possible to 
say which of the two coupling constants is larger. We remark that these two values were 
estimated by assuming G = 8.2 K, however, taking into account the uncertainty of ±1 K 
in the estimate of Curie- Weiss temperature and that exact diagonalizations provide useful 
estimates for J2/ Ji < 0.4, while QMC simulations only for J2/ Ji ^ 2I, it is difficult to assign 
an error bar to these estimates of J2/Ji- 

A more accurate determination of the ratio J2/J1 can be done by analyzing C"^{T) 
data in the light of diagonalization results by Singh and NarayananS and of the numerical 
calculations by Bacci et al.0. From the numerical results reported in Refs. 15 and 16 it is 
possible to plot the amplitude of the specific heat at the maximum C""(T^) as a function of 
the ratio J2/J1 (Fig. 14a). It turns out that the value C™(T^) = (0.436 ±0.004)i? found for 
Li2VOSi04 {R = NAks) (see Fig. 2) is compatible only with J2/J1 = 0.44±0.01 or l.liO.l 
(Fig. 14a). To discriminate among the two ratios one can analyze how = 3.5 ± 0.1 K 
varies as a function of J2/ Ji- Since 

^ = -f(l + 5^), (7) 

with T^/Q = 0.42 ± 0.04, one can check which value of J2/ Ji is compatible with the results 
of T^/Ji vs. J2/ Ji reported by Singh and Narayananlil (see Fig. 14b). One observes that 
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Eq. 7 is satisfied only for J2/J1 around 0.1 or 1.1. Tlierefore, tlie only solution which 
is compatible with the experimental values both of and C(T^) is J2/ Ji = 1-1 ± 0.1. 
This also indicates that G is close to 9 K (see Fig. 14b) and that — 0.59. Now, by 

assuming this value for one would derive from the analysis of the susceptiblity a value 

J2/J1 ~ 2, in agreement with the specific heat analysis, even if an accurate estimate with 
QMC is prevented, since in this range of J2/ Ji the results start to suffer from the minus sign 
problem. A value of J2/J1 around 1.1 also implies that Li2VOSi04 lies on the right hand 
side of the phase diagram reported in Fig. la, where the ground state is expected to be a 
collinear phase, in complete agreement with NMR results below Tc (see later on). 

Further information on the superexchange constants can be achieved from the analysis 
of '^Li 1/Ti. In the limit T ^ Ji + J2, 1/Ti is constant (see Fig. 12) and, by resorting to 
the usual gaussian form for the decay of the spin correlation function, one can writeS 

(l/TOco = ^^^^ X J: 141^ (8) 

with y4*^- (i,j = x,y,z) the components of the hyperfine tensor due to the k*^ V^"*" and 7 



the nuclear gyromagnetic ratio, uje = y Jl + J2{^b/K)^2zS{S + 1) /?> is the Heisenberg 
exchange frequency, where 2; = 4 is the number of nearest neighbour spins of a V^^ coupled 
either through Ji or J2. By using in Eq. 8 (1/Ti)oo = 0.2 ms~^, i.e. ''Li relaxation rate in 



the T range 300 — 3.2 K, one finds y + ~ 8.7 K, close to what one would derive from 
susceptibility and specific heat measurements. 

On decreasing temperature ^Li and ^^Si 1/Ti remain constant down to 3.2 K, at variance 
with yU+ 1/Ti (usually called A), which diverges on decreasing T, already at 4.2 K, due to the 
growth of the AF correlations. Both for nuclei and the spin-lattice relaxation is induced 
by the fluctuations of the effective local field, driven by the correlated spin dynamics, and 
one can write 

= A = ^ 5: \A^\^SUl ^r) (9) 

where l^g-p is the hyperfine form factor and Saa{<l,^^R) (« = x,y,z) are the components of 
the dynamical structure factor at the NMR or /iSR resonance frequency. One immediately 
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realizes that a different trend of NMR and /iSR 1/Ti can originate from the different form 
factors, which couple each one of these probes in a different way with the spin excitations at 
the critical wave- vector. In particular, one might suspect that ^Li and ^^Si form factors filter 
out the AF correlated spin excitations. However, if one considers that ^Li is coupled via a 
transferred hyperfine interaction with V^^ nearest neighbours (see Eq. 6), one finds that ''Li 
form factor is little g-dependent and that no filtering of the AF excitations can be envisaged. 
Moreover, the divergence of ^Li and ^^Si 1/Ti at Tc evidences that the fluctuations at the 
critical wave vectors cannot be completely filtered out. 

Another relevant difference is still present between /iSR and NMR measurements. While 
the former were performed in zero field, NMR 1/Ti measurements were carried out in mag- 
netic fields ranging from 1.8 to 9 Tesla, at which the Zeeman energy is comparable to the 
superexchange couplings. Therefore, it is tempting to associate the different behaviour of 
NMR and /iSR spin-lattice relaxation rates above Tc to a crossover of regime induced by the 
magnetic field. In particular, the T independent Ti measured in NMR would be consistent 
with a quantum critical regime^, while the exponential divergence of 1/Ti (A) measured 
with yuSR would be consistent with a renormalized classical regimeEl, where, by resorting to 
classical scaling arguments for 2D systems, one can write0 

T „ 



l/Ti(T) = A(T) oc ^(T) = 0.493a x 6^""=/^ 



1 -0.43j + O(j)2 



(10) 



with C, the in-plane correlation length, a the lattice step and ps the spin stiffness. From the 
T dependence of A (see Fig. 7) above Tc one derives 2ttPs = 7.4 K, less than the value 1.150 
expected for a non-frustrated systemB. 



B. Below Tc 

Since V^^ magnetic moments lie in the ab plane, as suggested by susceptibility measure- 
ments (Fig. 3) and by the EPR analysis of the g tensor^, and provided that the dipolar 
magnetic field cancels at ^^Si site, one realizes that the order must be collinear with a critical 
wave vector Q = (vr/a, 0), where x is the direction of the magnetic moment Jl. 
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The second order transition to the low T coUinear phase is evidenced by the peaks in 
1/Ti and in dx/dT. It is remarkable to observe that Tc is practically field independent up 
to at least 9 Tesla (see Fig. 4), where gfisH/kB > Ji + J2, while a decrease is expected, 
with Tc vanishing for gf^sH ~ 6A;b</2, if J2/ Ji — 1, i-e. at ~ 20 Teslal^. A possible 
explanation for this peculiar behaviour is that the structural distortion occuring just above 
Tc, deduced from ^^Si NMR spectra, causes an increase in the coupling constants and that 
even at 9 Tesla gfisH/kB < Ji + J2. Another possibility is that Li2VOSi04 is a 2D XY 
system with Tc close to the corresponding Berezinskii-Kosterlitz-Thouless transition^. 

Also the T dependence of the sublattice magnetization, derived either from the local field 
at the muon or from the splitting of ^Li NMR line, was found independent on the magnetic 
field intensity from zero up to 9 Tesla. From ZF/iSR measurements it has been possible to 
derive a critical exponent (3 = 0.235 ± 0.009 for the sublattice magnetization (see Fig. 6). 
Remarkably, this value of (3 is very close to the one predicted for a 2D XY model on a finite 
sizeS. Although some in-plane anisotropy can be discerned from the susceptibility data just 
above (see the inset to Fig. 3), there is no evidence of a crossover from Heisenberg to 
XY in the T dependence of the correlation length, derived from A(T) above tJhI. It is also 
interesting to observe that the sublattice magnetization measured by means of /iSR shows 
a slight high T tail, as expected in a finite size systemii. If the order is purely 2D, without 
long range order along the c axis, one would expect ^Li nuclei, which lie between V^"*" layers, 
to be characterized by a broad powder-like NMR spectrum. This is certainly not the case for 
Tc — T ^ 0.2 K (see Fig. 9), however one cannot exclude from the NMR measurements that 
the order is 2D in the very vicinity of T^. In fact, since the strong in-plane XY correlations 
enhance the 3D coupling the difference between the 2DXY and 3D ordering temperatures is 
expected to be small, of the order of the interlayer coupling J^B. An upper limit for J± can 
be estimated by assuming that Tc is the 3D ordering temperature of a Heisenberg AF, where 
Tc ~ OAJ±^^(Tc^. From the temperature dependence of /i"*" relaxation rate A (see Eq. 10) 
one finds C,(Tc)/a ~ 5.3, leading to J_l ~ 0.2 K. This value is possibly overestimated and 
difficult to justify if one considers the chemical bondings in Li2VOSi04 structure. Therefore, 
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a purely 2D order should be observable only for Tc — T ^ 0.2 K. 

Although the nature of this phase transition remains to be clarified, one can argue that 
the insensitivity both of the Neel temperature and of the critical exponent of the sublattice 
magnetization to the magnetic field indicates that the phase transition is driven by the XY 
anisotropy. 

In 3D magnets with two or more possible pitch vectors Q, the ordering usually cor- 
responds to a choice of pitch vector. The situation is often more complicated at lower 
temperature, and further transitions corresponding to other combinations of the pitch vec- 
tors or to the appearance of higher harmonics have been reported. Besides, the relevant 
parameter for the nature of the transition is the product N = n x m, where n is the number 
of components of the order parameter (3 for Heisenberg) and m is the number of equivalent 
wave-vectors0. As a consequence, the resulting transition can have a large critical exponent 
j3, typically around 0.4, or might in some cases be discontinuous. 

The results reported in the present paper suggest that the transition is split into two 
transitions: First a structural transition, as revealed by Si NMR, then an ordering transition, 
as seen at the Li site. A natural question arises as to whether the Ising degree of freedom 
corresponding to the two possible collinear states is associated with the structural distortion 
or with the magnetic ordering. We believe that the first possibility is the most likely both 
on experimental and theoretical grounds. Experimentally, the small value of the exponent 
/3 is typical of layered magnets with XY symmetry. If the parameter was increased by a 
factor 2 with respect to the number of components of the order parameter due to the Ising 
degree of freedom, one would not expect to observe such a small exponent. Besides, the 
choice of a pitch vector for the collinear phase renders the two directions inequivalent, and 
this is likely to be coupled to the lattice and to be associated with a structural distortion. 

One has to notice that the structural distortion occuring just above Tc may have modified 
the spin hamiltonian. Therefore, a discussion of the properties of the ordered phase on the 
basis of the parameters extracted above could be misleading. Nevertheless, one has to 
notice that, to be consistent with a collinear order, J2/J1 must be larger than ~ 0.65 also 
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below Tc- 

Information on the coupling constants below can be derived from the T dependence 
of ^Li nuclear spin-lattice relaxation. Below ^Li 1/Ti is mainly driven by two-magnon 
Raman processes^, leading to a T dependence for T ^ A, the gap in the spin-wave 
spectrum, and to 1/Tl oc T^exp(— A/T) for T ^ A. The low T dependence of ^Li 1/Ti 
turns out to be activated and, by fitting the data for T < 2.2 K with the latter expression, 
one finds A = 6 ± 1 kI. This value of the spin-wave gap is quite large if compared to 
the value of 6 = Ji -|- J2, estimated from susceptibility measurements above Tc, and would 
imply an axial anisotropy D ^ = 8.2 ± 1 K ~ A^/(Ji + J2)), which is quite large for 
V^^. In fact, the values of the g factor estimated from ESR measurements are very close to 
2 and yield a value oi D < \ kEiI. Moreover, if D ~ O Li2VOSi04 should behave as an Ising 
system, not as an XY or Heisenberg one, in sharp contrast to the experimental findings. 
Thus, one is tempted to argue that the low T collinear phase is characterized by coupling 
constants slightly larger than the ones determined above T^, so that D <ti Ji + J2 and its 
absolute value is smaller. 

Finally, one has to expect that frustration also causes a reduction of the staggered mag- 
netization due to the enhancement of quantum fluctuations. The T — average magnetic 
moment of V^"*" ions can be obtained from ^Li NMR spectra below Tc. By extrapolatimg to 
T ^ the splitting of ''Li NMR satellites and taking into account the hyperfine couplings 
given by Eq. 6, one can estimate a V^^ magnetic moment /i(T — > 0) ^ 0.24/^^. This value is 
reduced not only with respect to the value 0.65/1^ expected for a non-frustrated 2DQIIAF, 
but also with respect to the value derived numerically by Schulz et al.0 for J2/ Ji — 1, 
suggesting that probably below Tc J2/J1 ^ 1- 

IV. CONCLUSION 

In conclusion, it has been shown that Li2VOSi04 is a prototype of a frustrated 2DQHAF 
on a squre lattice with J2/J1 — 1.1 and J2 + Ji = 8.2 ± 1 K. Its ground state is a collinear 
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phase, as expected for J2/J1 ^ 0.65. The phase diagram as a function of the magnetic field 
intensity is characterized by a constant Tc{H), for < if < 9 Tesla. This observation, 
together with the fact that the critical exponent of the magnetization /3 ~ 0.235, suggest 
that the transition to the coUinear phase is driven by the XY anisotropy. The structural 
distortion occurring just above Tc, is expected to lift the degeneracy between the two coUinear 
ground states and to modify the superexchange couplings. In order to gain further insights 
on the nature of the phase transition and on the effective coupling constants below Tc further 
measurements with other techniques (e.g. inelastic neutron scattering) are required. 
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FIGURES 

FIG. 1. a) Schematic phase diagram of a frustrated 2DQHAF on a square lattice as a function 
of the ratio J2/J1 of the superexchange couphngs. b) Structure of Li2VOSi04 projected along 
[001]. Si04 tetrahedra are in gray, VO5 pyramids are in black while the gray circles indicate Li"*" 
position. For details see Ref. 9. 

FIG. 2. a) Temperature dependence of Li2VOSi04 molar specific heat below 70 K. The solid 
lines shows the phonon contribution to C(T), according to Eq. 1 in the text, with Qd = 280 
K. b) Magnetic contribution to the specific heat, obtained after subtracting the phonon term 
corresponding to the solid line in a) 

FIG. 3. Temperature dependence of the susceptibility % = M/H, for = 3 kGauss, in 
Li2VOSi04 powders. The dashed line shows the best fit according to Eq. 2, for 15 < T < 300 K. 
In the upper inset the derivative dx/dT is reported, evidencing a phase transition around 2.8 K. 
In the lower inset magnetization measurements in a Li2VOSi04 single crystal, both for H parallel 
and perpendicular to the c axis are reported. The intensity of M for ^ _L c have been rescaled for 
the sake of comparison. 

FIG. 4. Magnetic field versus T phase diagram for Li2VOSi04 . The circles indicate the field 
dependence of Tc derived from the kink in the susceptibility and/or from the peak in dx/dT (see 
Fig. 3), while the squares the corresponding values of Tc determined from ^Li NMR spectra (see 
Fig. 10). 

FIG. 5. Time evolution of /i+ polarization in Li2VOSi04 powders for T close to T^. The solid 
regular line shows the best fit according to Eq. 3 in the text. The T stability was within ±5 x 10~^ 
K. 

FIG. 6. Temperature dependence of the local field at the muon in Li2VOSi04 powders, derived 
from ZF/xSR measurements. The dashed line indicates the critical behaviour for a critical exponent 
of the magnetization (3 = 0.235 ± 0.009 (see text). 
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FIG. 7. Temperature dependence of the muon longitudinal relaxation rate in Li2VOSi04 pow- 
ders. The solid line indicates the T dependence of A according to Eq. 10, with a spin stiffness 
Ps = 7.4/27r K. 

FIG. 8. Plot of "^Li NMR paramagnetic shift versus the macroscopic susceptibility in Li2VOSi04 
, for H II c. The solid line shows the best fit yielding a total hyperfine coupling of 2.6 kGauss and 
a chemical shift (5 = —70 it 30 ppm, for ^ || c. 

FIG. 9. ^Li NMR spectra for = 1.8 Tesla along the c axis in a Li2VOSi04 single crystal, in 
the proximity of T(.. 

FIG. 10. Temperature dependence of the splitting of ^Li NMR satellites, for H = 1.8 Tesla 
along the c axis. The solid line shows the critical behaviour for an exponent (3 = 0.24. 

FIG. 11. 2^Si NMR powder spectrum in Li2VOSi04 for H = 1.8 Tesla in the proximity of Tc. 
The dotted lines mark the position of the peak at high and at low T. 

FIG. 12. ^Li nuclear spin-lattice relaxation rate l/Ti for H \\ c in Li2VOSi04 , for = 1.8 
Tesla (open squares) and 9 Tesla (closed circles). The dotted line gives the best fit according to 
the expression for 2-magnon relaxation processes (see text), yielding A = 6 it 1 K. In the inset the 
corresponding T dependence in the range 1.6 to 100 K is reported. 

FIG. 13. 29si NMR l/Ti in Li2VOSi04 for H = 1.8 Tesla, for T < 4.2 K. 

FIG. 14. a) Amplitude of the maximum in the molar specific heat for a frustrated 2DQHAF 
versus Jil J\- The open squares represent the data derived from Ref. 15, while the closed circles 
derived from Ref. 16. The gray region around C"^{T^)/R = 0.436 represents the experimental 
value for this quantity, inclusive of the error bar. b) T^/Ji (see text) versus J2/J1 derived from 
Ref. 15. The solid lines show the behaviour according to Eq. 7, for values of Q corresponding to the 
lower and upper limits of the Curie- Weiss temperature estimated from susceptibility measurements. 
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Fig.l (R. Melzi et al.) 
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